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ABSTRACT 


The influence of ionospheric refraction on orbit determination has been studied through 
the use of the Orbit Determination Error Analysis System (ODEAS). This paper presents 
the results of a study of the orbital state estimate errors due to the ionospheric refraction 
corrections, particularly for measurements involving spacecraft-to-spacecraft tracking 
links. 6 

In current operational practice at the Goddard Space Flight Center (GSFC) Flight Dy- 
namics Facility (FDF) , the ionospheric refraction effects on the tracking measurements 
are modeled in the Goddard Trajectory Determination System (GTDS) using the Bent 
ionospheric model. While GTDS has the capability of incorporating the ionospheric re- 
fraction effects for measurements involving ground-to-spacecraft tracking links, such as 
those generated by the Ground Spaceflight Tracking and Data Network (GSTDN), it does 
not have the capability to incorporate the refraction effects for spacecraft-to-spacecraft 
tracking links for measurements generated by the Tracking and Data Relay Satellite Sys- 
tem (TDRSS). The lack of this particular capability in GTDS raised some concern about 
the achievable accuracy of the estimated orbit for certain classes of spacecraft missions 
that require high-precision orbits. Using an enhanced research version of GTDS, some 
efforts have already been made to assess the importance of the spacecraft-to-spacecraft 
ionospheric refraction corrections in an orbit determination process. While these studies 
were performed using simulated data or real tracking data in definitive orbit determina- 
tion modes, the study results presented in this paper were obtained by means of covari- 
ance analysis simulating the weighted least-squares method used in orbit determination. 


The current operational version of ODEAS has the capability to compute ionospheric 
refraction corrections for range and range-rate measurements for both GSTDN and 
TDRSS data using the Bent ionospheric model. Using this capability in ODEAS, this 
study demonstrates how the magnitude and the characteristics of the spacecraft-to- 
spacecraft ionospheric refraction depend on such factors as spacecraft altitudes, solar 
activity, tracking geometry, and the local solar times of measurements. This study also 
provides error analysis results showing how the ionospheric refraction affects the space- 
craft state estimate errors for routine orbit determination of spacecraft missions such as 
Gamma Ray Observatory (GRO), depending on the tracking geometries and the tracking 
pass lengths. ODEAS is a software system particularly convenient for analyses of this 
type. The results obtained from this study will provide a quick reference for ionospheric 
refraction effects on orbit determination and will be useful for assessing orbit accuracy 
requirements for many future spacecraft missions. 


* supported by the National Aeronautics and Space Administration (NASA) /Goddard Space Flight Center 

iGSFCj, Greenbeit, Maryland, under Contract NAS 5 - 31500 . 
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1.0 INTRODUCTION 


The Earth’s ionosphere extends from approximately 80 kilometers to beyond 1000 kilometers from the sur- 
face of the Earth and consists of ionized particles. As a tracking signal travels through this medium, the 
refraction or bending of the signal path occurs due to the presence of free electrons. The magnitude of 
refraction is dependent upon the electron densities along the path of the tracking signal and the frequency of 
the signal. To achieve high accuracy in determining a spacecraft orbit, it is therefore necessary to properly 
correct the tracking data for the refraction effects. 

In the current operational practice at the Goddard Space Flight Center (GSFC) Flight Dynamics Facility 
(FDF) , the ionospheric refraction effect is computed in the Goddard Trajectory Determination System 
(GTDS) for ground-to-spacecraft tracking links using a global electron density model known as the Bent mode 
(References 1 and 2) . However, no refraction correction is performed for spacecraft-to-spacecraft inks for 
Tracking and Data Relay Satellite System (TDRSS) data; therefore, TDRSS tracking data are only partly 
corrected for ionospheric refraction. To circumvent this deficiency, a geometrical editing criterion can be 
employed in GTDS to reject the tracking data subject to significant ionospheric disturbance. This practice 
reduces the amount of tracking data actually used for the orbit determination process. It is therefore desirable 
to know the characteristics and the magnitude of the ionospheric refraction corrections and how they influence 
the accuracy of orbit determination. Such information will be useful to an orbit analyst in assessing the accu- 
racy requirements of the spacecraft missions and will provide insights into ways of achieving more accurate 
spacecraft orbits by minimizing the orbit errors associated with the ionospheric refraction corrections. Using an 
enhanced research version of GTDS, some efforts have already been made to assess the importance of the 
spacecraft-to-spacecraft ionospheric refraction corrections in an orbit determination process (Reference 3). 
While these studies were performed using simulated data or real tracking data in definitive orbit determination 
modes, the study results presented in this paper were obtained by means of covariance analysis simulating the 
weighted least-squares method used in orbit determination. 

This paper presents results obtained by using the Orbit Determination Error Analysis System (ODEAS). 
ODEAS is a general-purpose linear orbit determination error analysis system which can be used as a covariance 
analysis tool to analyze the orbit determination errors resulting from various systematic error sources and 
random noises. Throughout the study, the Earth’s ionosphere was modeled as a nonhomogeneous medium 
based upon the Bent model. The study results can be divided into two parts. The first part presents the 
magnitude and the general characteristics of ionospheric refraction correction as a function of spacecraft-to- 
spacecraft tracking geometry, solar activity, orbit inclination, user spacecraft orbital height, and local times of 
TDRS and the user spacecraft at the time of measurement. The second part investigates the spacecraft orbit 
determination errors that can result from neglecting the ionospheric refraction effects in spacecraft-to-space- 
craft links. Due to the linearity assumptions involved in the orbit estimation process, the results obtained can 
also be extended to address the effects of inexact modeling of ionospheric refraction in spacecraft-to-spacecraft 
tracking links by simple linear scaling. The results for a spacecraft mission such as the Gamma Ray Observa- 
tory (GRO) are discussed. 

A brief background of the Bent model is provided in Section 1.1 of this paper. Section 1.2 describes the 
mathematical procedures used in ODEAS for modeling the ionospheric refraction correction for spacecraft-to- 
spacecraft links using the Bent model. Section 2 discusses the study results, and Section 3 presents conclu- 
sions and recommendations. 

1.1 BENT MODEL 

The Bent model is an empirical world-wide electron density model named after its original developers 
R. B. Bent et al. (Reference 2). This model was derived in the early 1970s from analysis of the many 
thousands of satellite and ground-based ionospheric soundings and satellite electron density measurements 
available at the time. It was extensively evaluated before incorporation into GTDS (Reference 4). The mode 
provides the electron density profiles (electron density as a function of altitude) derived by incorporating the 
effects of variations in latitude, longitude, solar activity, geomagnetic activity, seasonal variations, and loca 
time. Using this model, an electron density profile similar to that shown in Figure 1 can be constructed for any 
geographic location around the globe. This profile is represented in the Bent model using seven curve sections, 
including two parabolic sections and five exponential tail sections. 
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ELECTRON DENSITY, N 


Figure 1. Empirical Worldwide Electron Density Profile 

The computation of the electron density profile consists of evaluating parameters specifying the upper and 
lower parts of the parabolic section including the maximum density, N m , the corresponding altitude, h m , and 
r^o P f° nenUal constants specifying the five exponential tail sections. This model has been successfully used in 
GTDS for correcting the ground-to-spacecraft tracking links for ionospheric refraction. The detailed mathe- 
matical description of the model as implemented in GTDS is described in Reference 1. Although the model 
has been successful in correcting the ground-to-spacecraft tracking measurements, the validity of the model for 
correcting spacecraft-to-spacecraft measurements, especially for high-altitude regions, has yet to be thoroughly 
evaluated for application in GTDS. 

The following section discusses the ionospheric refraction corrections using the Bent model for the spacecraft- 

to-spacecraft links as implemented in ODEAS. Reference 5 provides the mathematical details of the model 
implementation. 

1.2 IONOSPHERIC REFRACTION CORRECTION FOR SPACECRAFT-TO- 
SPACECRAFT LINKS 

The ionospheric refraction effects can be characterized in terms of the variable local index of refraction, n, of 
the medium through which the signal is propagated. This index of refraction can be expressed as 

n = 1 - N, (!) 


153 




where Nj is the refractivity and is given in terms of the electron density, N (electrons per cubic meter), and 
the signal transmission frequency, f (hertz), as follows: 

N, = 40.3 j (2) 

For a given measurement, the signal frequency, f, is a constant, and the electron density, N, varies with the 
altitude, longitude, latitude, and local time of the signal path. The range correction due to ionospheric refrac- 
tion, Api , can be computed by integrating Ni as expressed in Equation (2) along the signal path, as follows: 


40.3 

a 6i = £2 


I 


N ds 


(3) 


where ds denotes an infinitesimal ray path length. 

For a spacecraft-to-spacecraft tracking link, the integral expressed in Equation (3) must be evaluated by a 
numerical method, because the electron density, N, cannot easily be expressed as a function of the ray path, s. 
In ODE AS, this integral is evaluated using a 12-point Gaussian quadrature scheme. 

Figure 2 illustrates the two tracking configurations for a spacecraft-to-spacecraft tracking link that are relevant 
for the effects of ionospheric refraction. For configuration 1, (ft > 90 degrees), the integral of Equation (3) 
is evaluated from the user spacecraft at U to the maximum ionospheric height at Q, which is set at 3000 kilo- 
meters. For configuration 2, (£ <90 degrees), the integral is evaluated in two parts divided at point P, the 
point of lowest approach for the signal path. One integration is done from point P to Q and the other from 
point P to U. 




(a) CONFIGURATION 1 <p > 90°) 

Figure 2. Spacecraft-to-Spacecraft Segment (Two Possible Configurations Shown) 


In the ODEAS implementation, to achieve computational efficiency, a new electron density profile may not 
always be recomputed at each integration point along the ray path. A new electron density profile is computed 
for point Sj whenever the geographic longitude and latitude of this point are different from those of the 
previous integration point on the ray path, Sj_i, by more than a specified tolerance, i.e., 

I* - *,-il + |Ai - Am I £ Tol (4) 
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where 


geographic latitude and longitude pairs for the points s* and Sj_i , 
respectively 

user-specified tolerance level 


(0 i( Ai) and (0 M , Am) = 


Tol 


Therefore, by adjusting the tolerance level, the number of electron density profiles computed in the integration 
can be controlled. Throughout the study, the tolerance level for recomputing the electron density profile was 
set to 1 degree. This setting ensures recomputation of electron density profiles at each integration point along 
the ray path for most of the tracking geometries except for those where the integration points are very close to 

each other, such as in the case of large f} angles, in which case the range ionospheric corrections obtained by 
setting Tol = 1 degree and Tol = 0.001 degree differed by less than 2 percent. 

Range-rate corrections can be computed from range corrections as follows: 


A q _ Agi(tj) - Apifti-!) 

' At 


where 

Agi(tj) and Ap^tj.]) = two successive range corrections computed at time tj andt^ , respectively 

At = time difference between two measurements (usually the Doppler count 

interval) 


2.0 RESULTS 

The results of this study are presented in two parts. Section 2.1 discusses the general characteristics of iono- 
spheric refraction correction as a function of various parameters. Section 2.2 describes the orbit determina- 
tion error analysis results that can be expected by ignoring the ionospheric refraction effects in 
spacecraft-to-spacecraft tracking links using the routine orbit determination scenario of a GRO-type spacecraft 
mission as an example, 

2.1 GENERAL CHARACTERISTICS OF IONOSPHERIC REFRACTION CORRECTIONS 

As stated earlier, the magnitude of ionospheric refraction correction is a function of tracking signal frequency 
and the electron densities along the signal path as it traverses the ionosphere. This is expressed mathematically 
in Equation (3). For a given measurement, the frequency, f, is a constant; and the ionospheric refraction 
correction is complicated due to the nature of electron density representation. This section presents the 
magnitude and the general characteristics of ionospheric refraction correction as a function of various parame- 
ters that affect the electron densities. The parameters studied were limited to spacecraft-to-spacecraft tracking 
geometry, solar activity, orbital height and inclination, and local solar times of TDRS and the user spacecraft at 
the time of measurement. Although the ionospheric refraction can also be affected by the seasonal variation, 
this effect was not studied parametrically. 

To study the effects of these parameters, 18 different specially designed circular orbits were investigated using 
ODEAS. These basic spacecraft orbits were constructed by the combination of six different orbital heights 
(350, 550, 750, 950, 1150, and 1350 kilometers) and three different inclinations (28.5, 63.14, and 99.03 de- 
grees). tor each basic run scenario, a 24-hour definitive period was simulated using TDRSS two-way range 
and range-rate tracking. The tracking signals were modeled as S-band signals, with a frequency of 2100 mega- 
hertz, and a data sampling rate of 30 seconds. As shown in Equation (3), the ionospheric refraction correc- 
tion is proportional to the inverse of the signal frequency squared; and, therefore, the results obtained can be 
scaled appropriately for different signal frequencies. For example, the ionospheric refraction correction value 
obtained with a K-bana signal would yield a value approximately 50 times smaller than the S-band signal. 
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A measurement was scheduled whenever the user spacecraft was visible to the TDRS. This was done to ensure 
the collection of a fairly large sample of tracking data with a variety of tracking geometries. Ionospheric 
refraction corrections were computed for all the spacecraft-to-spacecraft range and range-rate data. The re- 
sults obtained are described below. 

2.1.1 Tracking Geometry 

The spacecraft-to-spacecraft tracking geometry can be characterized by the relative positions of the relay 
spacecraft (TDRS), the user spacecraft, and the Earth. In the TDRSS environment, the relay spacecraft is in 
geosynchronous orbit at a height of approximately 35,800 kilometers above the Earth’s surface. Figure 2 
illustrates the two tracking geometries. In Figure 2a, both the relay and the user spacecraft are on the same 
side of the Earth; whereas in Figure 2b, the relay and the user spacecraft are on the opposite sides of the 
Earth. 

The elevation angle of the relay with respect to the user, a , can be defined as the line-of-sight elevation angle 
of the relay spacecraft as measured from the local horizontal plane at the user. This elevation angle can be 
computed by 


a - p - 90 degrees 


( 6 ) 


With this definition, the elevation angle, a , assumes a positive value in Figure 2a and a negative value in 
Figure 2b. 

The height of ray path, h 0 , is the height of a point of closest approach on the ray path. For Figure 2a, h 0 is 
the same as the user spacecraft orbital height. The central angle, <5 , is defined as the angle between the 
position vectors of TDRS and the user spacecraft, subtended at the center of the Earth. 

It can be seen in Figure 2 that as the elevation angle, cl , decreases, there is a greater possibility that the ray 
path will traverse a more “diverse” ionosphere, which consists of many different electron density profiles due 
to bigger possibility of variations in latitude, longitude, and local times. Conversely, as the elevation angle 
increases, the ray path encounters the ionosphere with less diversity. Figures 3a and 3b illustrate the average 
ionospheric refraction for range and range-rate data as a function of elevation angle for various orbital heights 
at a 28.5-degree inclination. The ionospheric refraction correction value plotted at a given orbital height at a 
given elevation angle was obtained as an average value in the case of range correction, or as a root-mean- 
square (RMS) value in the case of range-rate correction, from the collection of tracking data taken at that 
specified elevation angle during a 24-hour period. RMS values were used for range-rate data because the 
ionospheric refraction corrections for range-rate data could assume positive and negative values, and the RMS 
can represent the mean magnitude of the correction regardless of sign. 

Figures 4a, 4b and Figures 5a, 5b illustrate the same type of plots for orbital inclinations of 63.14 and 
99.03 degrees. To see the spread of ionospheric refraction correction values obtained at each elevation angle, 
Figures 6a and 6b illustrate the maximum, minimum, and average or RMS values of the ionospheric correc- 
tions for a 28.5-degree inclination, 350-kilometer altitude orbit. 

All these figures demonstrate the fact that the ionospheric refraction correction for both range and range-rate 
measurements decreases as the elevation angle increases. The elevation angle dependence is more pro- 
nounced for lower orbital heights (350 and 550 kilometers) and less pronounced for higher altitudes 
(1150 and 1350 kilometers). 

Figures 3 through 5 also illustrate the functional dependence of ionospheric refraction corrections on the 
orbital height. It can be seen from these figures that the magnitude of corrections are generally higher for 
lower orbital heights. At a 350-kilometer altitude, the average ionospheric refraction corrections for range 
measurements are about 1.5 to 2 times as high as those values observed at a 550-kilometer altitude. As the 
orbital height increases, the correction values become smaller and its dependence on the orbital height become 
less significant. For the spacecraft orbital heights of 950, 1150, and 1350 kilometers, the corrections are 
found to be in the same neighborhood. 
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Daily Solar Flux : 218.9 
12-Month Average 

Solar Flux : 200.6 
Sun Spot Number : 163 

Signal Frequency : 2100 MHz 
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Elevation Angle (Degrees) 

Figure 3a. Average Range Ionospheric Correction Versus Elevation Angle for 
a 28.5-Degree Inclination for Various Altitudes 


Daily Solar Flux : 218.9 
12-Month Average 

Solar Flux : 200.6 
Sun Spot Number : 163 

Signal Frequency : 2100 MHz 
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Figure 3b. RMS Range-Rate ionospheric Correction Versus Elevation Angle for 
a 28.5-Degree Inclination for Various Altitudes 
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Daily Solar Flux : 218.9 
12-Month Average 

Solar Flux : 200.6 
Sun Spot Number : 163 

Signal Frequency : 2100 MHz 


Elevation Angle (Degrees) 

Figure 4a. Average Range Ionospheric Correction Versus Elevation Angle for 
a 63.14-Degree Inclination for Various Altitudes 
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0-5 0 5 10 15 20 25 30 3b 

Elevation Angle (Degrees) 

Figure 4b. RMS Range-Rate Ionospheric Correction Versus Elevation Angle for 
a 63.14-Degree Inclination for Various Altitudes 


158 










Average Range Ionospheric Correction (Meters) 


9 - 



Elevation Angle (Degrees) 

Figure 5a. Average Range Ionospheric Correction Versus Elevation Angle for 
a 99.03-Degree Inclination for Various Altitudes 



Elevation Angle (Degrees) 

Figure 5b. RMS Range-Rate Ionospheric Correction Versus Elevation Angle for 
a 99.03-Degree Inclination for Various Altitudes 
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Figure 6a. Minimum, Maximum, and Average Range Ionospheric Correction Versus Elevation 
Angle for a 28.5-Degree Inclination and a 350-Kilometer Altitude 



Figure 6b. Minimum, Maximum, and RMS Range-Rate Ionospheric Correction Versus Elevation 
Angle for a 28.5-Degree Inclination and a 350-Kilometer Altitude 
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frhi?r?r Cd Tit^’ t ^ e *^ PPer b °“ nd ° f the ionos P here was set at 3000 kilometers. This choice is somewhat 
nr^ntln • A ?J 0U * h h ® ^ ontnbutions from altitudes above 3000 kilometers are not included in the results 

way t^the TDRS Sud’e LT® 'T* ““ * eStimated by evaluatin 8 the integral of Equation (3) all the 
way to the TDRS altitude. Analyzing these magnitudes obtained for a few sample tracking passes it was found 

that they could vary significantly according to the altitudes of the user spacecraft. For aTacecraft at 
a i ude of 350 kilometers, for example, the contribution from above 3000 kilometers to the measurement 
correction was, in general, less than 5 percent of the contribution from below 3000 kilometers. However for 
a spacecraft a 1350 kilometers, the contribution from above 3000 kilometers was found to be aTwgh a 
40 percent of the contribution from below 3000 kilometers. 8 " 

2.1.2 Local Solar Times of Measurement 


Hmt i0 T°hu he ff iC denSity f ° r any location over the Earth ’ s surface is also a function of the local solar 

time. This effect has been incorporated in the Bent electron density model. Because of this functional de- 
pendence, it is conceivable that the ionospheric refraction corrections for a particular range Grange-rate 
measurement can be dependent on the local solar times of the relay and user spacecraft. For the spacecraft 
orbits studied, the ionospheric refraction correction was found to be affected more significantly by die local 
times of user spacecraft, with maximum corrections occurring around local times of 12 to 16 hours. This effect 
was consistently observed for all the spacecraft orbits studied. 


Figure 7 shows the average range ionospheric refraction corrections as a function of the local solar time of the 

ZSSTtu J h V Pa “ Craft orbiK •» « <■«. differenTorWtal kiclina- 

' r "l 28 ' 5, 63 ‘! 4, and 99 03 degrees ~ at a fuced alti tude of 350 kilometers. This figure indicates a big in- 
se in ionospheric refraction correction values around the user spacecraft local times of 12 to 16 hours 
peaking at approximately 15 hours. The same trend was also observed for RMS range-rate corrections 



“^ t Pbe ^ me "° n Rested by the Bent model closely matches the real ionospheric behavior, this may 
uggest that to effectively exclude the tracking data subject to high ionospheric disturbance, it might also be 
necessary to take into account the local time of the user spacecraft as a data editing criterion, together with the 
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tracking geometry considerations. A data editing criterion based on tracking geometry characterized by the 
height of ray path and the central angle is presently employed in the operational GTDS. In this practice 
tracking data obtained at a height of ray path below a certain prescribed value and a central angle greater than 
a certain prescribed value are excluded. The prescribed values used are specific to mission requirements. 

2.1.3 Solar Activity 

To study the effects of solar activity on the ionospheric refraction correction, three sets of solar data were used, 
as summarized in Table 1. FI, F2, and F3 represent cases of moderately high, medium, and low solar activi- 

ties, respectively. 


Table 1 . Solar Data Used for Ionospheric Refraction Simulations 


SIMULATION 

NUMBER 

DAILY 

SOLAR FLUX* 

12-MONTH AVERAGE 
SOLAR FLUX* 

SUN SPOT 
NUMBER 

FI 

218.9 

200.6 

163 

F2 

159.2 

111.4 

56 

F3 

74.4 

76.1 

20 


* UNITS = 10- 22 WATTS/METER S /HERTZ 


Figures 8 and 9 illustrate the average ionospheric refraction correction of range data as a function of elevation 
angle and solar activity for 350-kilometer and 1350-kilometer altitudes, respectively, at a 63 14-degree 
inclination. It can be seen that the ionospheric refraction correction values increase as the solar activity 
increases. For range corrections as shown in Figure 8, the values for the FI curve are about two “mes higher 
than the corresponding F2 curve, and F2 is about two times larger than the corresponding F3 curve. Similar 
trends were also observed for RMS range-rate corrections. 

To study the sensitivity of ionospheric refraction corrections to the daily and 12-month average or mean solar 
flux values, range corrections for two tracking passes, one for a 350-kilometer orbit and one for a 1350-kilo e 
ter orbit, were computed using different daily and mean solar flux values. The daily and mean solar flux values 
were varied in steps of 50 units. Whenever the mean solar flux was changed, the corresponding sunspot 
number was also changed according to an empirical formula that relates the sunspot number to the mean solar 

flux (Reference 3). 

The results are summarized in Table 2. The values in cases 1, 3, and 5 were obtained using a common daily 
solar flux level of 218.9 with different mean solar flux levels. The values in cases i 2, 3, and 4 were obt 
with the same mean solar flux level of 200.6 and different daily solar flux levels. These results indicate that, in 
general, the refraction correction depends more strongly on the mean solar flux than on the daily flux, 
general trend, however, does not appear to hold for the part of a pass with lower elevation angles. 

2.1.4 Orbital Inclination 

From Figures 3 through 5, it can be determined that the average or RMS values of ionospheric refractlon 
obtained at various elevation angles for low-inclination (28.5-degree) orbits tend to e ig er t an 
high-inclination (99.03-degree) orbits. This tendency is found to be more pronounced in the low-altitud 
region (350 kilometers). At the 1350-kilometer altitude, the ionospheric refraction corrections for all thre 
orbital inclinations are in the same neighborhood. 

The observed dependence of the measurement refraction corrections on orbital inclination could be, at least 
partially, due to the strong dependence of the maximum electron density (N m ) and the correspon mg ®'g 
(h m ) on the latitude. Table 3 presents the values of these two parameters as a function of latitude. Th|se 
results were obtained using the Bent model at a longitude of 318 degrees and a loca so ar ime o 
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Average Range Ionospheric Correction (Meters) 3 Average Range Ionospheric Correction (Meters) 



8. Average Range Ionospheric Correction Versus Elevation Angle for a 63.14-Degree 
Inclination, a 350-Kilometer Altitude, and Various Flux Levels 



Elevation Angle ([Degrees) 


Figure 9. Average Range Ionospheric Correction Versus Elevation Angle for a 63.14-Degree 
Inclination, a 1350-Kilometer Altitude, and Various Flux Levels 
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Table 2. Daily and Mean Solar Flux Dependence of Ionospheric Refraction 
Corrections for Range Measurements 


EPOCH: 

ORBIT: 


FEBRUARY 1, 1990; 0 h 0 m 0 s 
350-KILOMETER ALTITUDE; 28.5-DEGREE 


INCLINATION 


TIME 

f MINI! ITPS 

ELEVATION 

ANGLE 

RANGE CORRECTIONS (METERS) 

FROM EPOCH) 

(DEGREES) 

CASE 1 

CASE 2 

CASE 3 

CASE 4 

CASE 5 

5 

7.6 

15.15 

20.23 

24.05 

28.23 

24.42 

10 

24.6 

2.35 

3.64 

4.20 

4.87 

5.76 

15 

41.4 

1.08 

1.61 

1.79 

1.97 

2.52 

20 

55.6 

0.80 

1.14 

1.26 

1.39 

1.73 

25 

61.3 

0.74 

1.02 

1.13 

1.25 

1 .55 

30 

53.7 

1.22 

1.70 

1.96 

2.20 

2.77 

35 

38.8 

2.00 

2.92 

3.58 

4.34 

5.12 

40 

22.0 

1.58 

2.14 

2.47 I 

2.74 

3.11 

45 

5.0 

5.10 

7.04 

8.13 

9.35 

10.01 

DAILY SOLAR FLU 

X* 

218.9 

168.9 

218.9 

268.9 

218.9 

12-MONTH MEAN SOLAR FLUX* 

150.6 

200.6 

200.6 

200.6 

250.6 


EPOCH: FEBRUARY 1, 1990; 0 h 0 m 0‘ 
ORBIT: 1350-KILOMETER ALTITUDE; 


28.5-DEGREE INCLINATION 


TIME 

/MINUTES 

ELEVATION 

ANGLE 

RANGE CORRECTIONS (METERS) 

FROM EPOCH) 

(DEGREES) 

CASE 1 

CASE 2 

CASE 3 

CASE 4 

CASE 5 

5 

2.6 

0.61 

mm 

mm 

1.33 

* I 

10 

16.4 

0.18 

W&m 


0.34 


15 

30.1 

0.10 

0.13 

0.13 

0.14 

0.17 

20 

43.3 

0.07 

0.09 

0.09 

0.09 

0.11 

25 

54.6 

0.07 

0.08 

0.08 

0.08 

0.09 

30 

60.5 

0.07 

0.08 

0.08 

0.09 

0.09 

35 

57.4 

0.08 

0.09 

0.09 


0.11 

40 

47.4 

0.10 

0.12 

0.14 

0.17 

0.21 

45 

34.6 

0.12 

0.14 

0.15 

0.18 

0.19 

50 

20.9 

0.15 

0.18 

0.19 

0.19 

0.21 

DAILY SOLAR FLU 

X* 

218.9 

168.9 

218.9 

268.9 

218.9 

12-MONTH MEAN SOLAR FLUX* 

150.6 

200.6 

200.6 

200.6 

250.6 


•SOLAR FLUX UNITS = lO -22 WATTS/METER 2 /HERTZ 


As can be seen from this table, both h m and N m , especially the latter, increase in the lower latitude region. 
The parameter N m acts as a scale factor for the entire electron density profile, thus affecting the electron 
densities computed in all different height regions. The variation of the maximum electron density height ( m ) 
can also contribute to an increase in the refraction corrections for orbits with lower inclinations. The values'o 
N m and h m as a function of latitude were examined for several other longitudes, and all showed the trend 
similar to that described in Table 3. 
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Table 3. Latitude Variation of the Maximum Electron Density ( N m ) 
and the Corresponding Height (h m ) 


LATITUDE 

(DEGREES) 

I'm 

(KILOMETERS) 

N m 

(10’° ELECTRONS/ 
METER 3 ) 

LATITUDE 

(DEGREES) 

hm 

(KILOMETERS) 

N m 

(10'° ELECTRONS/ 
METER 3 ) 

88.47 

82.38 
73.69 
61.99 

51.39 
44.31 
34.30 
21.03 
19.38 

382.13 

363.65 

342.92 
325.73 
319.81 
320.52 
328.01 

352.92 
398.08 

60.708906 

49.976265 

65.968391 

118.187695 

170.900019 

192.604224 

201.218013 

213.455625 

272.916559 

-88.40 

-78.17 

-70.16 

-59.07 

-47.15 

-38.04 

-30.12 

-18.44 

-7.19 

379.49 

361.63 

352.95 

346.42 

345.39 

349.87 

361.66 

409.45 

464.41 

33.860825 

45.700007 

54.567174 

79.233256 

124.645481 

164.401967 

196.020009 

215.195575 

170.750405 


NOTES: EPOCH = FEBRUARY 1 , 1990 

DAILY SOLAR FLUX = 218.9 ( 10' 22 WATTS/METER 2 /HERTZ) 

12-MONTH MEAN SOLAR FLUX = 200.6 ( 10* 22 WATTS/METER 2 /HERTZ> 
LONGITUDE = 318 DEGREES 
LOCAL SOLAR TIME = 15 h 


2.2 


™« CTS ° F NEGLECTING SPACECRAFT-TO-SPACECRAFT IONOSPHERIC 
REFRACTION ON ORBIT DETERMINATION ACCURACY 


As stated in Section 1, the current operational practice using GTDS ignores the ionospheric refraction effects 
or the spacecraft-to-spacecraft links of TDRSS measurements. It is therefore desirable to know for different 
ypes of spacecraft missions the magnitude and characteristics of the spacecraft orbit errors that can result bv 

IhaMnAnenr! 5?*?' ^ *5* 10n ° Spheric refraction effect is functionally dependent upon various parameters 
t nfluence the electron density, a comprehensive study of this kind would require extensive investigations 

nvolvmg combinations of various mission scenarios and various parameters. This section presents the orbit 
ermmation error analysis results for the routine orbit determination scenario of a GRO-type spacecraft. 

Because the primary interest of the study was the ionospheric refraction effects, the study was limited to 
^° t W . ] , °" os P henc ^fraction of spacecraft-to-spacecraft links influences the orbit determination 
Reference 6 ^ ° rblt yS ‘ S lnvesti 8 ation for the GRO spacecraft mission using ODEAS can be found in 

2-2.1 Error Analysis Results 

° R ° spacecra ft is scheduled to be launched in late 1990. In this study, the nominal operational orbit of 
GRO spacecraft was modeled as a near-circular orbit, with an orbital height of 350 kilometers and an 

28-5 d ® grees ' . The orbital P eriod for this orbit was 91.54 minutes. Routine orbit determination 
r ® GR ° spacecraft was simulated by processing 30-hour definitive data arcs in which the orbital state of the 
u spacecraft and the atmospheric drag scaling factor ( p, ) were simultaneously estimated together with the 
orhi't 3 6 em ® nts TDRS-East and TDRS-West. To see the effects of ionospheric refraction on the spacecraft 
errors for different tracking geometries and tracking pass lengths, four tracking scenarios (T1 through T4) 

Tnpee , a u Cd: thC det3ilS ° f theSC scenarios are § iven in Table 4. Ail four tracking scenarios involved one 
K66 tracking pass per GRO revolution, alternating between TDRS-East and TDRS-West. For scenarios T1 
and T2, tracking pass lengths were limited to 5 minutes each; whereas in scenarios T3 and T4, tracking pass 
engths were limited to 20 minutes each. In scenarios T1 and T3, the tracking measurements are subject, to 

re atively high ionospheric refraction effects; whereas in scenarios T2 and T4, the ionospheric refraction effects 
are relatively small. 
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Table 4. Simulated Tracking Scenarios for GRO Routine Orbit Determination 
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of the error ana |ysis simulations, the ionospheric refraction corrections for TDRS-East-to-GRO and 
DRS-West-to-GRCJ tracking links were included as considered error sources, with their respective 3(7 errors 
S K, t0 b J 100 pe “; ent of spacecraft-to-spacecraft ionospheric refraction correction. Error analysis results 
obtained for the four tracking scenarios are summarized in Table 5. The contribution from the range and 

InTrnnc w eaS f Ure . mentS arC ^en separately for spacecraft-to-spacecraft tracking involving the TDRS-East 
rel ay spacecraft. Total contributions from the range and range-rate measurements listed in 
the table are obtained by treating these quantities as fully correlated parameters. The maximum position errors 
listed are the maximum errors encountered during the 30-hour data arc for each position component error 
which may not necessarily occur at the same time. V 

Error analysis results for scenarios T1 and T2 show that the GRO spacecraft position errors contributed by 
spacecraft-to-spacecraft ionospheric refraction can vary significantly, depending on the geometry of the track- 
ing measurements included m the orbit estimation process. As shown in Table 4, the number of tracking 
passes and measurements included in T1 and T2 are the same, except that the measurements included in T1 
have substantially lower elevation angles than those of T2 and are therefore subject to higher ionospheric 
disturbance than T2. Because of this, the GRO spacecraft position errors contributed by the spacecraft-to- 
spacecraft ionospheric refraction are found to be substantially higher in T1 than in T2, as shown in Table 5. 
T n n ” a . no T1, the maximum root-sum-square (RSS) position errors contributed by the TDRS-East and 
TDRS-West tracking measurements are 57.82 and 91.90 meters respectively; whereas in scenario T2, the same 
position error contributions are only 3.78 and 9.76 meters, respectively. The results suggest that by properly 
selecting the tracking measurements with favorable tracking geometries, the maximum RSS position errors can 
bejeduced by as much as 50 to 80 meters. The same trend can be observed by comparison of scenarios T3 


Comparison of scenarios T1 with T3 and T2 with T4 also show the effects of tracking pass lengths on the 
spacecraft position errors. In the tracking scenarios studied, the spacecraft-to-spacecraft ionospheric refraction 

it t CtS 6 GR ? ® pacecraft position errors are found to be generally smaller for longer tracking pass lengths, 
t s ould be noted that in all the simulated cases, the tracking scenarios associated with the longer tracking pass 
lengths also have better tracking geometries. Therefore, the effects observed in the above scenarios may be 
due to the combined effect of longer tracking pass lengths and better geometries. 


To illustrate the relationship between the ionospheric refraction correction versus the resulting spacecraft posi- 
tion errors, Table 6 lists the ionospheric refraction corrections and the corresponding maximum RSS position 
errors contributed separately from the range and range-rate measurements. A general trend can be observed 
rom this table indicating that an increase in the average value (for range measurement) or the RMS value (for 
range-rate measurement) of ionospheric refraction correction results in an increase in the maximum RSS posi- 
tion errors. However, this general trend does not hold very well between scenarios T4 and T2 where the 
magnitudes of corrections are small. 


3.0 CONCLUSIONS AND RECOMMENDATIONS 


The magnitudes of ionospheric refraction correction for spacecraft-to-spacecraft tracking links can vary de- 
pending upon the various parameters that affect the electron densities. The parameters studied include the 
local solar times of the measurements; solar activity; spacecraft inclination; and the tracking geometry, which is 
characterized by (1) the elevation angle of the TDRS with respect to the TDRS user spacecraft and (2) the 

orbital height. The following general conclusions can be made regarding the spacecraft-to-spacecraft iono- 
spheric refraction: 

* T ^ e ionospheric refraction corrections for both range and range-rate measurements increase with 
the decrease in the elevation angle of the TDRS with respect to the TDRS user spacecraft. This suggests that to 
avoid high ionospheric disturbance in tracking data, data taken at very low elevation angles should be ex- 
cluded. For high spacecraft altitudes (950 kilometers and up), the magnitudes of ionospheric corrections are 
much smaller than those for the lower altitudes, but the overall trend of elevation angle dependence still exists. 

• The ionospheric refraction correction is found to be a strong function of the local solar time of the 
TDRS user spacecraft at the time of measurement. It was found that measurements taken between 12 to 
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Table 5 GRO Spacecraft Position Errors Contributed by Spacecraft-to-Spacecraft Ionospheric 
Refraction Effects for Different Simulated Tracking Scenarios 



•i < 
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Table 6. Relationship Between the Spacecraft-to-Spacecraft Ionospheric Refraction 
Corrections Versus the Corresponding GRO Spacecraft Position Errors 


TRACKING* 

SCENARIO 

RELAY 

SPACECRAFT 

AVERAGE 

IONOSPHERIC 

CORRECTION 

RANGE 

(METERS) 

MAXIMUM 
POSITION 
ERROR 
DUE TO 
RANGE 
(METERS) 

RMS 

IONOSPHERIC 
CORRECTION 
RANGE RATE 
(MILLIMETERS/ 
SECOND) 

MAXIMUM 
POSITION 
ERROR 
DUE TO 
RANGE RATE 
(METERS) 

T1 

TDRS-EAST 

18.2 

18.25 

98.66 

60.09 

T3 

TDRS-EAST 

7.97 

5.01 

51.55 

21 .93 

T4 

TDRS-EAST 

3.05 

1.08 

3.95 

2.80 

T2 

TDRS-EAST 

2.84 

1.15 

2.59 

2.95 

T1 

TDRS-WEST 

17.1 

17.07 

84.99 

96.17 

T3 

TDRS-WEST 

8.57 

5.12 

48.20 

21.64 

T4 

TDRS-WEST 

4.70 

1.86 

10.08 

6.62 

T2 

TDRS-WEST 

4.14 

1.90 

7.15 

8.39 


T1 AND T2 = 5-MINUTE TRACKING PASS PER GRO REVOLUTION, ALTERNATING BETWEEN 
TDRS-EAST AND TDRS-WEST 

T3 AND T4 = 20-MINUTE TRACKING PASS PER GRO REVOLUTION, ALTERNATING BETWEEN 
TDRS-EAST AND TDRS-WEST 


16 h° urs of local solar time tend to have significantly higher ionospheric disturbance than those taken at 
different local solar times. If this phenomenon manifested by the Bent model closely matches the real iono- 
spheric behavior, this may suggest a need to take into account the local time of the user spacecraft as a data 
editing criterion, together with the tracking geometry considerations, to effectively exclude the tracking data 
subject to high ionospheric disturbance. A data editing criterion based on tracking geometry, characterized by 
the height of the ray path and the central angle, is presently employed in the operational GTDS. To verify this 
local solar time effect, future studies should include assessing the orbit determination accuracy using real track- 
ing data to determine whether more accurate spacecraft orbits can be realized by using a data editing criterion 
that excludes those tracking data taken at a user spacecraft local solar time around 15 hours and at very low 
elevation angles. 

• The ionospheric refraction correction is found to be a strong function of solar activity. For a 
spacecraft orbit at a 350-kilometer altitude and 28.5-degree inclination, the ionospheric refraction correction 
values obtained at a^ moderately high solar activity [average solar flux = 200.6, daily solar flux = 218.9 (units = 
10 watts/ meters /hertz)] can be twice as large as the corrections obtained at a relatively lower solar activity 
[average solar flux = 111.4, daily solar flux = 159.2 (units = 1(T 22 watts/ meters 2 /hertz) ] . It was also found 
that the ionospheric refraction effect is more sensitive to the 12-month average solar flux level than the daily 
solar flux level. 


* The maximum electron density value, N m , as defined on the Bent electron density profile, has a 
strong latitude dependence. It was found that N m values tend to be higher in low latitude (near-equatorial) 
regions than in the high-latitude (polar) regions. It is therefore conceivable that in a particular spacecraft orbit, 
t e ionospheric refraction effects may tend to be high if many of the tracking measurements were taken when 
the user spacecraft was near the equatorial region. This may explain why the ionospheric refraction corrections 
observed m the study tend to be generally higher for the low-inclination (28.5-degree) orbits than those of the 
high-inclination (99.03-degree) orbits. 


From the error analysis results of the GRO spacecraft, the following conclusions can be drawn: 


vary signi 


The spacecraft position errors contributed by spacecraft-to-spacecraft ionospheric refraction c 
nficantly, depending on the geometry of the tracking measurements included in the orbit estimation 
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process. In a routine orbit determination involving 30 hours of data with 5 minutes of TDRSS tracking per 
GRO revolution, selecting the tracking measurements with favorable tracking geometries can reduce the space- 
craft maximum RSS position errors from 57.82 meters to 3.78 meters (in the case of TDRS-East-to-GRO 
ionospheric refraction contribution) and from 91.90 meters to 9.76 meters (in the case of TDRS-West-to-GRO 
ionospheric refraction contribution). 

• An increase in the mean value of the ionospheric refraction correction is found to be associated 
with an increase in the maximum RSS position errors. However, this general trend does not hold very well 
where the magnitudes of the ionospheric refraction corrections are small or the geometric quality of the proc- 
essed data is appreciably altered. This general relationship also implies that parameters affecting the mean 
ionospheric refraction values, such as solar activity, will also affect the corresponding spacecraft orbit errors in 
a similar fashion. 
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